62 research outputs found
An unconditionally stable finite difference scheme for solving a 3D heat transport equation in a sub-microscale thin film
AbstractHeat transport at the microscale is of vital importance in microtechnology applications. The heat transport equation is different from the traditional heat diffusion equation since a second-order derivative of temperature with respect to time and a third-order mixed derivative of temperature with respect to space and time are introduced. In this study, we develop a finite difference scheme with two levels in time for the 3D heat transport equation in a sub-microscale thin film. It is shown by the discrete energy method that the scheme is unconditionally stable. The 3D implicit scheme is then solved by using a preconditioned Richardson iteration, so that only a tridiagonal linear system is solved for each iteration. The numerical procedure is employed to obtain the temperature rise in a gold sub-microscale thin film
Enhancing accuracy for solving American CEV model with high-order compact scheme and adaptive time stepping
In this research work, we propose a high-order time adapted scheme for
pricing a coupled system of fixed-free boundary constant elasticity of variance
(CEV) model on both equidistant and locally refined space-grid. The performance
of our method is substantially enhanced to improve irregularities in the model
which are both inherent and induced. Furthermore, the system of coupled PDEs is
strongly nonlinear and involves several time-dependent coefficients that
include the first-order derivative of the early exercise boundary. These
coefficients are approximated from a fourth-order analytical approximation
which is derived using a regularized square-root function. The semi-discrete
equation for the option value and delta sensitivity is obtained from a
non-uniform fourth-order compact finite difference scheme. Fifth-order 5(4)
Dormand-Prince time integration method is used to solve the coupled system of
discrete equations. Enhancing the performance of our proposed method with local
mesh refinement and adaptive strategies enables us to obtain highly accurate
solution with very coarse space grids, hence reducing computational runtime
substantially. We further verify the performance of our methodology as compared
with some of the well-known and better-performing existing methods
A new compact finite difference scheme for solving the complex Ginzburg-Landau equation
a b s t r a c t The complex Ginzburg-Landau equation is often encountered in physics and engineering applications, such as nonlinear transmission lines, solitons, and superconductivity. However, it remains a challenge to develop simple, stable and accurate finite difference schemes for solving the equation because of the nonlinear term. Most of the existing schemes are obtained based on the Crank-Nicolson method, which is fully implicit and must be solved iteratively for each time step. In this article, we present a fourth-order accurate iterative scheme, which leads to a tri-diagonal linear system in 1D cases. We prove that the present scheme is unconditionally stable. The scheme is then extended to 2D cases. Numerical errors and convergence rates of the solutions are tested by several examples
- …